Pair condensation of bosonic atoms induced by optical lattices 
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We design a model of correlated hopping for bosonic atoms in optical lattices. Such model 
exhibits three kinds of phases, comprehending a Mott insulator, a charge density wave and a pair 
quasi-condensate. The implementation of the model relies on two-state atoms embedded in state- 
dependent lattices and having spin-dependent interactions. Contrary to other models of pairing, 
our design is not based on perturbative effects and should be observable in current experiments. 



Since the achievement of Bose-Einstein condensation 
in alkali atoms [H, 0, 01 , we have witnessesed two major 
breakthroughs in the many-body physics of cold atoms. 
One is the realization of Cooper pairing and the BCS 
to BEC transition with fermions [3, H, B]- The other 
one is the implementation of lattice Hamiltonians using 
neutral atoms in optical lattices 0, S|- Supported by 
this success, many theoretical papers suggest using cold 
atoms with two goals: the quantum simulation of well 
known Hamiltonians such as Hubbard models § and spin 
lattices [§] , and the quest for new physics such as bosonic 
quantum Hall effect [3, U| and lattice gauge theories 
13 , 13j |. In this work we aim at the latter, introducing 
a robust mechanism of pairing that works for bosonic 
atoms, opens theoretical challenges and is suited for the 
recent experiments in optical superlattices [3, EH • 

Pairing is a central concept in strongly correlated 
states. In particular, it is the essence of ordinary BCS su- 
perconductivity. In BCS theory, an attractive interaction 
mediated by a phonon bath is the basis by which elec- 
trons merge into bosonic pairs that conduct electricity 
without friction. Similar phenomena appear in the BCS 
pairing of fermionic atoms, where the natural attraction 
is enhanced by Feschbach resonances [1, [H, 0] . 

Another, less known mechanism for pairing is corre- 
lated hopping. It appears naturally in fermionic tight- 



binding models [1 
turn magnetism 1 2 




and in quan- 



, consisting on the motion of particles 
being influenced by the environment. This is normally 
reflected by terms of the form nia^ak appearing in the 
Hamiltonian. Correlated hopping could lead to the for- 
mation of bound electron pairs 3, 13] and it has been 
put forward as an explanation for high T c superconduc- 
tivity [HEH. 



We will introduce a different mechanism for pairing, 
which is based on transport-inducing collisions. As 
shown in Fig. [1] when atoms collide they can mutate 
their internal state. If the atoms are placed in a state- 
dependent optical lattice, whenever such a collision hap- 
pens, the pair of atoms must jump to a different site 
associated to their new state. For deep enough lattices, 
as in the Mott insulator experiments [7| , this coordinated 
jump of pairs of particles will be the dominant process 



and the atoms will become a superfluid of pairs. 

This work consists of three parts. The first one begins 
with a setup of cold atoms in a deep optical superlattice. 
We show how collisions combined with state-dependent 
lattices lead to a strong correlated tunnelling, as de- 
scribed above. The outcome is an effective Hamiltonian 
which can be related to the experimental parameters. 
We will then draw a realistic phase diagram using both 
Gutzwiller [H and Matrix Product States (MPS) [27| 
variational methods. Both procedures essentially iden- 
tify three kinds of phases: a charge-density wave (CDW), 
a Mott insulator (MI) and a pair condensate or pair su- 
perfluid (PSF). In the last part of the paper we consider 
different experimental procedures to detect the different 
insulating and superfluid orders. 

Let us begin by considering the setup in Fig. [Ik, where 




FIG. 1: (a) Two atoms in state |+) collide and change into 
state | — ). (b) When trapped in a state-dependent optical lat- 
tice, changing the state implies also jumping between lattice 
sites. We plot some forbidden (dashed) and allowed (solid) 
processes. Only pairs of atoms can hop between lattices, (c) 
In order to maximize the overlap between atoms in different 
states and thus the strength of correlated hopping, in this 
paper we consider two state-dependent superlattices. 



bosonic atoms in two internal states, a+ and a_, are 
trapped in state-dependent superlattices that spatially 
overlap. The hopping outside a double well is negligible. 
In this limit, we have two bands of localized states with 
energies proportional to the hopping amplitude on a dou- 
ble well, +J and — J. In particular the lowest band will 
be given by the bosonic states c 2 k '■= -^{o-2k+ + o-2k+i+), 

and c 2 fc+i := -^(a 2 k+i- + a%k+2-)- 

If the interaction is weak compared to the band gap 
2J, only the coupling between states in the same band 
will be relevant. We will focus on a particular situation 
in which the trapped states are actually dressed states, 
a± := (af ± a|)/\/2, and the interaction is diagonal 
in those other operators, a-\ and oj,. This situation is 
reached by combining the optical lattice configuration in 
Ref. [28| with asymmetric on-site interactions 



U 



: (n T +n l ) 2 : + U^n^ 



(1) 



Note that the term Uin^n^ is not invariant under rota- 
tions of the internal state of the atoms. When expressed 
in the |±) basis, it is this term that scatters pairs of atoms 
from one cloud to another [Fig. [T]a,]. Since the superlat- 
tices have a relative displacement, these collisions give 
rise to transport. Using the approximation Uq, U± <C 2J 
one arrives at our central model 



H 



U : (m + tij) : + V riiUj 
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t c] Cj 



(2) 



where U = (2U Q + U 1 )/16, V = -Ux/8, and t = Ui/16. 

Out of these terms, V favors neighboring sites to be 
filled, t is the correlated hopping that spreads pairs of 
atoms through the lattice, and U > restricts the total 
number of atoms preventing collapse. If 17 = and V = 
0, we expect the bosons to form a superfluid of paired 
particles, while for strong interactions there should be a 
quantum phase transition to an insulator. 

The first evidence of pairing is obtained by studying a 
lattice with only two particles. This setup has two po- 
tential ground states: one with both particles isolated 
and zero energy, and a delocalized state of a single pair 
Impair) oc £ 4 cf |0) W ith energy E 2 = 2(U - V). There- 
fore, if U < V or Uq < Ui/2, the pair of bosons will be 
bound, E 2 < 0. An important question is whether this 
criterion survives in the presence of other atoms, or if 
many-body effects may destroy the pairing. 

In order to study the many-body physics of our model 
we have used the Gutzwiller ansatz, which is a variational 
method based on the product state [26[ 
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Minimizing the expectation value of the free energy, 
F := H — fiN, with respect to the variables /„ under 



the constrain of fixed norm, \f n | 2 = 1, one obtains 
the phase diagram in the phase space of interaction and 
chemical potential, (Ui,U ,[i). For the sake of simplicity 
we have used U\ := 1 as unit of energy everywhere. 

The results are shown in Fig. [21 where we plot the 
expectation values of the density, n = (n) = (e'e), the 
particle number variance, An := y (n 2 — n 2 ), and what 
we identify as the order parameter of the paired super- 



fluid, 



From the zeros of the density fluctuations we 



recognize the Mott regions with n = 1, 2 and 3 particles 
per site. The boundary of the region with n = 1 can 
be estimated analytically using a Gutzwiller state with 
nonzero components /o, /i and f 2 , which gives 



U 



, and fi h = U 



(4) 



These are the energies to add a particle or to make a 
hole in this insulating lobe. If fi < /i/j, the Gutzwiller 
ansatz gives fractional densities but, as we will see below, 
this is an artifact of the uniform trial wavefunction. If 
fj, > j-ip and Uq < Ui/2, we obtain a region of nonzero 
(c 2 ). We take this as a sign of a long-range coherence in 
the two-body density matrix which cannot be attributed 
to an ordinary Bose-Einstein condensate, since the same 
simulation gives (c) = 

The Gutzwiller ansatz does not reproduce accurately 
neither the location of phase transitions nor the behav- 
ior of correlators. For instance, C\ := (c\ci + A) and 
:= ( c 1 2c 1+a) computed with Eq. ([3]) are both uni- 
form functions. In order to study these properties, we 
have used the MPS method (27J to estimate variationally 
the ground state of our model. The MPS is a more com- 
plex wavefunction that, for large enough computational 
resources and not too large systems, should reproduce 
the correlations in the superfluid regime. Calculations 
with up to 40 sites reveal that matrices of size D = 40 
are enough to pinpoint the different phases. 

In Fig. [3] we plot the same quantities as for the 
Gutzwiller ansatz. The contour plots reveal that the 
n = 1 and n = 2 insulating regions still exist, with a 
similar size and shape as in the Gutzwiller case. As we 
explain later, the region above /i p again shows a paired 
superfluid behavior, which extends even beyond the limit 
Uq = Ui/2 where the binding energy of an isolated pair 
becomes positive [Fig. [3b] . The biggest difference is in 
the triangle < fx < [Ah- I n this region of low density, 
the bosons arrange forming a charge density wave, i. e. 
a pattern alternating and 1 atom per site. 

Regarding the superfluid, quasi-long-range order is 
identified by a slow decay of the off-diagonal elements 
of the single- and two-particle density matrix. As shown 
in Fig.[U the single particle correlator C\ — 5a, o is only 
different from zero at A = 0, where it becomes the den- 
sity. This could have indicated the presence of a Mott 
phase, were it not for the nonzero value of An (Fig. [SJs) 
and of the two particle correlator, C A , that decays slowly 



at long distances. Due to the size of our simulations, we 
have not yet been able to determine the behavior of C\ , 
but numerical fits of curves like Fig. |4] suggest an alge- 
braic decay, C\ ~ A~ a , with a nonuniversal exponent 
around a < 0.7 that depends on Uq and fi. 

Let us now introduce another possible experimental 
setup, illustrated in Fig.rjL and where two lattices trap- 
ping the two different species are shifted by half a pe- 
riod or A/4. As usual, the trapped states and the in- 
teraction are described in different basis. One can ex- 
pand the interacting Hamiltonian in terms of Wannier 
functions together with operators defining the trapped 
states. This leads to Eq. ([2]) with: U = (2U + Ui)/4, 
V = -(20U O + U x )/2 and t = -Ui/4. Here, U a and U x 
are the on-site interactions in the atomic basis and the 
factor p = J \w(x)\ 2 \w(x - X/4)\ 2 dx/ J {w^dx < f, 
computed using the Wannier functions w(x), measures 
the relative strength of interaction between different lat- 
tices with respect to those on the same lattice site. 

The phases that appear in our setup can all be recog- 
nized experimentally. First of all, the MI and the CDW 
have both a well defined number of particles per site and 
no coherence. Their time of flight pictures will have no 
interference fringes and the noise correlation will ex- 
hibit peaks at certain momenta 2^, 33] • In addition, the 
CDW corresponds to a setup where either n + or n_ are 
uniformly zero. Finally, the energy gap in these insula- 
tors can be probed by static 0] or spectroscopic means 

Regarding the pair superfluid, it is a perfect "conduc- 
tor" with a gapless excitation spectrum. Lacking single- 
particle order, C\ ~ a.e., it will also not show inter- 
ference patterns in the time of flight images. In order 
to measure C\ and detect the pairing, we suggest to 
use photoassociation and to build molecules out of pairs 
of atoms. Since the molecules will be built on-site, the 
nonzero correlator C\ will translate into long-range order 
for the molecules, which should now exhibit an interfer- 
ence pattern in time-of-flight images, slightly blurred due 
to the phase fluctuations inherent to ID. 

Finally, all energy scales arc referenced to the asymme- 
try of the interactions, U\. Therefore, the pair binding 
energy can be potentially as strong as the strength of the 
on-site interaction in a typical Mott insulator [7| , about 
1 kHz. In practice, while the most commonly used atoms 
such as Rubidium have smaller asymmetries, they can be 
enhanced using Feschbach resonances. 

Summing up, we have introduced a mechanism by 
which bosonic atoms with repulsive interactions can ex- 
hibit correlated hopping and pairing. The model @ ex- 
hibits multiple phases, among which the most relevant is 
a superfluid of paired bosons. All phases are connected 
by second order quantum phase transitions and can be 
produced and identified using variations of current exper- 
These ideas can be generalized to higher 
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FIG. 2: Ground state properties estimated with the 
Gutzwiller wavefunction (J3J). Grayscale plots of (a) the den- 
sity, n, (b) density fluctuations, An, and (c) pair condensate 
order parameter (c 2 ). Dashed lines mark the analytical esti- 
mates ((4]), while solid lines delimit regions of integer filling. 
All plots cover the same region (Uo, fi) £ [0, 5t/i] x [0, 10C/i] . 




FIG. 3: Ground state properties estimated with the MPS 
method. We plot (a) the fluctuations of the density, An, and 
(b) the averaged two-particle correlator, ^ A \C\\/L, where 
C\ := (cj 2 c 2 +A ). Solid lines delimit regions of integer filling. 
All plots cover the same region (Uo,n) € [0, 2Ui] x [0, 2Ui]. 
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FIG. 4: Single particle, C\ := (c^Ci+A) (dashed), and two 
particle, C\ := (a-l 2 af +A ), correlator (solid) for fx/Ui = 
0.9, Uq/Ui = and a lattice with 40 sites. The inset repeats 
the plot for C\ in log-log scale, with a fit C\ oc A ' 32 (dash). 



